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Abstract 

The program was ccmducted to determine the relative importance of entropy 
and direct caad>ustion noise in Lurbopropulsion systems and to determine the 
parameters upon idiich these noise sources depend. Theory and experiment were 
employed to determine that, at least with the apparatus used here, entropy 
noise can dominate combustion noise if there is a sufficient pressure gradient 
terminating the combustor. Measurements included combustor interior fluctuat- 
ing pressure, near and far field fluctuating pressure, and combustor exit 
plane fluctuating temperatures, as well as mean pressures and temperatures. 
Analysis techniques included spectral, cross-correlation, cross power spectra, 
and ordinary and partial coherence analysis. Also conducted were combustor 
liner modification experiments to investigate the origin of the frequency 
content of combustion noise. Techniques were developed to extract non-propa- 
gational pseudo-sound and the heat release fluctuation spectra from the data. 
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Introduction 


The purpose o£ this program was to Investigate the fundamental causes of 
core noise In turbopropulslon systems. The program was Initiated as M/ISA Con- 
tract NAS3-17861 for which Ref, (1) was Issued as a final report. The past 
three years of effort have been carried out under NASA Grant No, NSG 3015, 
and this report stunmarlzes those efforts. 

It was known at the program outset that there were at least two probable 

( 2 ) 

causes for core noise - entropy or Indirect noise and direct combustion noise. 
There had been no prior program aimed at Isolating the two suspected causes 
and determining their relative importance to the core noise problem. Moreover, 
there were some fundamental difficulties In understanding of the frequency 
content of direct combustion noise, as discussed in Ref, (1). This report con- 
tains experimental and theoretical Information aimed at separation of the ef- 
fect of entropy and direct combustion noise and resolution of the cause of the 
frequency content of direct combustion noise. 

The issue of the importance of core noise to the aircraft noise problem Is 
not addressed here. It Is clear that core noise represents a noise floor In 
current turbopropulslon systems but there is controversy concerning the strength 
of core noise relative to other sources^,^^*^^^* ^^^It is sufficient to remark here 
that core noise exists and Is measurable. This report will also add to current 
controversy because a major finding has been that entropy noise Is dominant 
In the apparatus used here, under conditions simulating those that would be 
found In engine installed configurations. Such a noise source has been rejected 
in Ref. (4) and correlation of engine noise results based upon a direct combus- 
tion noise theory has been claimed. While the experimental results reported 
here may be apparatus dependent, the theory suggests that both entropy and com- 
bustlon noise must be contended with In actual engines; moreover, the two sources 
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are currently linked and, hence, are not independent. 

This report is written in sunmary form with the details contained in 
appendices. Ibe conclusions contain reconiaepdations for further work, because 
other noise sources were found in the course of the work which may be important 
to the core noise problm. 
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Sunmarv 

The basic combustor unit used in the experimental work was taken from a 
Boeing 502-7D gas turbine unit. The unit is described in Ref. (1) and in 
Appendix A, It is a can type combustor with head end and side slots and is 
operated by discharging directly to the atmosphere (1 atm tests) or through 
a nozzle to the atmosphere (< 2 atm tests). In Appendix A the operating 
characteristics. such as efficiency and split between bypass and core flow, 
are discussed. Two can liners, other than the production liner, were manu- 
factured to vary the mean hydraulic diameter of the liner holes over a range 
of 4:1. These liners are also described in Appendix A. Hie hole size variation 
was introduced to test the hypothesis that the integral (macro) scale of the 
turbulence should play a role in the frequency of combustion generated noise. 

Types of measurements made on the combustor unit were far field noise, 
near field noise, interior pressure fluctuations, and exit plane temperature 
fluctuations and mean temperature. The types of analysis applied to the re- 
corded signals were spectral, cross spectral, cross correlation, ordinary co- 
herence and partial coherence analysis^^^ 

A theory of the noise behavior of combustors is developed in Appendices 
B and C. The theory is based upon the linearized vortf city-acoustic field 
equationi^ind several noise sources are recognized for the current experimental 
situation. Two of these which were not desired in the program were jet noise 
from the exhaust and the second was turbulence (or hydrodynamic or scrubbing) 
noise measured by the interior pressure transducer. These sources were not 
treated in the theory but were experimentally investigated. Another noise 
source which was eliminated from the theory was lateral vorticity - nozzle 

/Q \ 

interaction noise,'' ' This source, caused by velocity fluctuations entering 
a nozzle pressure gradient, must be investigated in a future program because 
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of the findings of the current program. However, study of the source was beyond 
the scope of the current program. The major emphasis in the theory was on direct 
combustion noise, a dilatation of the flow caused by fluctuations in the aggre- 
gate heat release, and entropy noise, pressure waves generated when hot (or 
cold) spots are convected through a nozzle pressure gradient. 

Entropy noise (as well as nozzle-vorticity noise) is absent if the combustor 
is exhausted directly to the atmosphere. Hence, the direct discharge tests were 
direct condsustion noise studies. Appendix B is concerned with these tests and 
the investigation of hydrodynamic noise contamin8.tion of interior pressure meas- 
urements. The primary cause of entropy fluctuations is temperature fluctuations 
entering a region of pressure gradient, such as a nozzle. The theory Indicates 
that for the apparatus used (^ 12% rms temperature fluctuations entering the 
nozzle and correlation length scales of the order of a centimeter) one should 
expect entropy noise to be seen experimentally. Ihis is especially true if a 
choked nozzle is attached to the combustor, as it was in the experiments of 
Appendix D. 

The tecbT'iques of ordinary and partial coherence analysis were used to make 
the determination of the relative dominance of combustion of entropy noise. 

Both techniques show the conclusion of entropy noise dominance at sufficiently 
high nozzle contraction. Partial coherence was, however, to have been the more 
sensitive Indicator of signal contamination by some other cause, such as vor- 
tlcity-nozzle interaction noise. Unfortunately, partial coherence analysis is 
highly sensitive to errors in experimental determination of spatial correlation 
of tenperature fluctuations. Since these correlations were not obtained with 
high accuracy, the accuracy of the partial coherence analysis suffered. The 
result was inconclusive determination of the existence of non-considered noise 
sources. 
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!Qie connection between entropy noise and combustion noise is investigated 
in Appendix E, It is shown there, both theoretically and experimentally, that 
the two noise sources can be either directly correlated or uncorrelated, de- 
pending upon the frequency range. For frequencies above the frequency defined 
by the inverse of a fluid particle stay time the two noise sources are in- 
coherent* Below this frequency there is a high correlation between them. The 
root cause, of course, is that both noise sources <.re related to the same heat 
release fluctuations. There is also evidence presented in Appendix D that a 
resonance feedback oscillation tal-es place at a frequency equal to the Inverse 
of a particle stay time. This is a feedback oscillation which links pressure 
fluctuations to entropy generation. While it takes place at a unique frequency 
and is of fundamental interest, it is of little interest in noise generation 
because of the small amount of acoustic power generated. 

Of interest from the data reduction standpoint is a novel method ^Ich was 
developed for measurement of thermocouple response time situ . This method is 
described in Appendix F, 

Substantial early work on this program was expended in development of an en- 
tropy pulse facility to measure the entropy admittance coefficient of various 
nozzles under various flow conditions. Results are only now forthcoming from 
this work and will be reported later in the archive literature. 

Finally, Appendix 6 addresses the issue of the frequency content of combus- 
tion noise. Ihe liner hole size variations were found to have virtually no effect 
upon the frequency content of the noise when exhausting directly to the atmos- 
phere. Both raw spectra of the interior microphone and derived spectra of the 
heat release fluctuations were used in the determination. This result was sur- 
prising, since it had been suspected that the mean length scale of turbulence, 
which should have been varied in these tests, was a strong factor in the fre- 
quency content.^ ^ Ihe frequency content issue bears further examination. 
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Conclusions 


It oust be borne In mind that some of the experimental conclusions here 
may be dependent upon the apparatus and conditions Investigated, Most Important 
among these were a) operation at 2 atm or below, b) use of a combustor with 
temperature fluctuations of about 12% of the mean absolute temperature and c) 
length scales of the temperature fluctuations of about 1 cm. 

1. Theoretically, using experimental Input, entropy noise will overtake and 
dominate combustion geucrated noise as a) the contraction ratio of a 
terminating nozzle Is Increased at fixed mass flow and fuel/alr ratio 

or b) the mass flow at fixed fuel/alr ratio Is Increased for a fixed 
contraction ratio. This takeover and dominance should be seen first at 
the higher frequencies (in the range investigated). 

2. Experimentally, conclusion # 1 Is borne out. 

3. At frequencies below roughly 200 Hz, entropy and combustion noise are 
highly correlated, since they have the same physical origin. 

4. Both Interior and exterior pressure measurements are contaminated by 
one or more noise sources, other than entropy or combustion noise. It 
was Inconclusive whether or not vortlclty-nozzle Interaction noise was 

' a dominant contaminant, and further work Is needed to Investigate this 
source. 

5. Measurement of the time-space cross power spectra of the temperature 
fluctuations must be performed with high accuracy If used In partial 
coherence analysis. In future work It Is recommended that a net of 
many thermocouples be used Instead of the two thermocouple traverses 
used here. 

6. Hydrod 3 rmunlc noise, ^Ich Is non propagatlonal pseudo-sound, dominates 
flush mounted Interior microphones below about 150 Hz. This may be 
largely eliminated by use of the infinite tube method with the micro- 
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phone displaced from the flow field. 

7. The frequency content of cond>u8tion noise is virtually unaltered by a change 
in the hydraulic diameter of the combustor liner holes. Past programs have 
found no combustor or flow variable that significantly alters the frequency 
content; however, fuel reactivity has not been varied over a wide range and 
deserves investigation in a future program. 
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Appendix A 

COIBUSTOR CHARACTERIZATION STUDIES 

The combustor used for die present noise measurement program was removed 
from a Boeing 502-7D gas turbine engine. This unit was designed to perform 
efficiently at the three-atmosphere, 41 o\ discharge of a single stage centrl> 
fugal compressor. The current burning apparatus supplies air at near ambient 
pressure and temperature to the cosdnistor. In addition, fuel nozzle flows are, 
at most, 40% of specified design values. Since off-design operation can signif i- 
cantly reduce burning efficiency, a performance evaluaticm of die combustor was 
judged necessary. As a first step, the air and fuel metering systems were checked 

3 

for accuracy. In the present blow-down facility, air is supplied from a 28 M 
storage vessel initially at 690 kPa (gage). Flow rates are determined at a regu- 
lated pressure with an orifice meter using mercury manometers to read upstream 
pressure and An orifice coefficient was computed from Ref. (Al) for the 
range of flows to be employed in the proposed noise test program. A check on 
computed air flows was made later in Che program during cold-flow pitot surveys 
at the combustor eidiaust plane. Agreement between Che two methods was found to 
be good. 

Fuel flows are determined with a Potter turbine meter and readout device. 

The unit was calibrated in water flow tests and was found to yield a linear 
output over its useful range. Maximum measured flow agreed within 0.5% of the 
manufacturer's stated value. The combustor unit was then mounted on the water 
flow test stand for fuel nozzle evaluation. The spray pattern was observed to 
consist of a barely discemable central plume, roughly 2 Qh in maximum diameter 
by 2.50m long, surrounded by an evenly distributed cone of fog-sized droplets. 
The outer spray cone intersect<id the liner wall about midway between the inlet 
ana exhaust planes. No significant alteration In spray pattern was seen over the 
range of flows used in Che current test program. In all cases, the nozzle spray 
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obseivatloas uere made with no air flow. 

A series of cold flow tests were next completed to determine the portion 
of air which completely bypasses the liner and eidiausts throng the annulus 
(Figure Al). Pitot surveys of the exit plane at the locations shown in Figure 
A1 for flows of 14.2 and 19.8 M^/Min. revealed that 25% + 1% of the total air 
mass flow was bypassed. For both cases the siun of the pitot-measured liner 
and annulus flows agrees with the orifice-measured flow within + 2%. Typical 
dynamic pressure profiles are shown in Figure A2 for the 14.2 M^/MLn. case. Pro- 
files for 19.8 M^/Miiu exhibit the expected higher pressure magnitudes but 
closely resemble those for 14.2 M^/Min. in distribution indicating no signifi- 
cant alteration for changing flow rate. 

Two cases with combustion reflecting a change in overall fuel/air ratio 
were selected for initial study; 14.2 M^/Min. air, (F/A)q^ = 0.016; 19,8 M^/Min, 
air, (F/A)q^ = 0.008. As in the cold flow tests, pitct measurements were per- 
formed at the stations shown in Figure Al. In addition, temperature measurements 
were made with a chrome 1 -a lumel thermocouple at the same locations. Mass flow 
calculations from these measurements using mean values of temperature and dy- 
namic pressure showed that tne by-pass flow was not altered appreciably when 
combustion was present. By-pass for both cases was 26% + of the total in- 
itial air flow. Total mass flow compmted frcrni the sum of measured liner and 
annulus flows in both cases was about 2,5% lower than that calculated from 
orifice measurements. This slight discrepancy probably reflects the diffi- 
culties in obtaining accurate dynamic pressure measurements in a reacting 
flow. Figures A2 and A3 present dynamic pressure and temperature profiles 
for the 14.2 M^/Mln. case with combustion. As in the cold flow tests, the 
dynamic pressure profiles with combustion for 19,8 M^/Min. are quite similar 
in shape to those for 14,2 M^/Min. and differ only in magnitude. The previous 
statement is also true when temperature profiles for the two cases are compared. 
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Figure Al, Combustor Test Apparatus 
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It should be noted that the temperature profiles reflect the earlier described 
fuel spray pattern In which a central plume Is enclosed In a conical skirt. 

Burner efficiencies were computed for the 14.2 M^/Mln« and 19.8 M^/Mln. cases 
with combustion. In each case, efficiency was calculated both with a thermody- 
namlc heat balance and with an expected temperature rise based on the adiabatic 
flame temperature. Each method Ignored the 26% annular flow In accordance with 
test results showing only minor heating. Figure A4 presents the computed effl> 
clency values and shows that the combustor performance was poor. It was observed 
from these results, however, that efficiency significantly improved with Increased 



overfall fuel to air ratio. On this basis. It was decided to Investigate combustor 
performance with Increasing fuel flow for a given selected air flow. 

Figure A5 shows the increase in measured combustor centerline exit temperatures 
for Increased fuel rates at an air flow of 14.2 M^/Min. Lines of constant effici- 
ency based on a thermodyitamic heat balance are also plotted for comparison. It 
should be noted that the mean exhaust temperature rise for these tests Is about 
80% of the centerline temperature rise so that comparison of experimental data 
with the efficiency curves is useful for determining trends only. Figure A5 
clearly Indicates improved performance as fuel flow Increases with diminishing 
returns for (F/A) 0 ^ values greater than 0.016. From these results It was decided 
to select an over-all fuel-air ratio of 0.021 for a more detailed exit plane 
velocity-temperature scan. 

The dynamic pressure and temperature profiles for (F/A)q^ « 0.021 are shown 
in Figs. A2 and A3. Figure A2 Indicates that d}mamlc pressure values decreased 
significantly from the (F/A)q^ * 0.016 case. It was observed during the (F/A)q^* 
0.021 tests tliat the flame routinely approaches and intermittently extends beyond 
the combustor can exit, indicating that conbustion Is not yet complete at the 
plane of measurement. As suggested earlier, measured dynamic pressure magnitudes 
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in reacting gases are unreliable, especially In this particularly severe case. 
Hie resulting calculated liner mass flows were, as expected, considerably low. 

A second test was performed to check the Initial measurements. Extra care was 
taken In obtaining the annular data since liner measurements were not considered 
reliable. Results of the second test verified the initial findings. By-pass flow 
was found to be about 227. of the Initial total air flow for this case. This 
small change from previous tests may Indicate an actual alteration of the com- 
bustor flow pattern, but probably reflects the physical difficulties Involved 
In making measurements In the annular flow with a flame extending past the exit 
plane. As predicted by the centerline temperature measurements of Fig. A5, the 
efficiency values for this condition were significantly Improved. The perform- 
ance Is tabulated In Fig. A4. 

The test program described above reveals that the present combustor unit 
must be operated at over-all fuel-air ratios on the order of 0.016 or greater 
to approach design efficiency. Accordingly, all subsequent measurements were 
performed at "high efficiency" flow ratios so that conclusions about combustor 
acoustics would not be subject to question from off-deslgn operation. 

Two alternate combustion liners were manufactured to study the effects of 
primary and secondary bleed hole width on combustor noise generation. In each 
unit axial slot area distribution and slot length was held equal to that of 
the original liner. With this constraint, machining and symmetry considerations 
dictated a minimum slot width equal to one-half that of the original ("S.w.x%") 
and a maximum slot width equal to twice that of the original ("S.W. x 2"). Area 
distribution was maintained by doubling the number of narrow slots and halving 
the number of wide slots (Including primary air Inlet slots surrounding the fuel 
nozzle on the liner dome). Hie new units were fabricated from duplicate liners 
obtained from the manufacturer. Each liner was modified by filling existing 



slots where necessary and milling new slots where required, 13ie original and 
alternate liners are shown schematically in Fig. A6, Bleed hold area distri- 
bution is plotted £or the three liners in Fig, A7. 

Operating characteristics of both liners were checked prior to the in- 
itiation of noise tests. It was observed that the S,W. x2 liner behaved much 
like the original liner in starting and steady state operation. Combustion 
appeared to be less complete than with the original liner as evidenced by 
larger areas of yellow flame at the liner exit plane. Flame distribution 
also was observed not to be as uniform as with the original liner. This ob- 
servation was confirmed in a post-test investigation when areas of discolora- 
tion on the liner surface where found. A check of exit plane temperatures at 
five locations for (F/A)q^ = 0,016 (air flow = 14.2 M^/Min) showed only a 
sli^tly lower mean temperature than for the original liner. 

Exploratory tests with the S.W. x% liner revealed that starting or main- 
taining a flame was impossible for over-all fuel to air ratios less than 0,0145 
at an air flow of 14,2 M^/Min, In contrast, the original and S,Wx2 liners 
could routinely operate at (F/A)q^ » 0.008 and lower. Improved distribution of 
air flow in the primary combustion zone for this liner apparently dilutes the 
fuel-air mixture uniformly below the flamability limits for the lower fuel 
flows. Unlike the original and S,w«x2 liners, the S,u.x% liner showed no 
evidence of yellow flame in the combustion zone. Ihe light-blue flame was well 
distributed at the exit plane and no discoloration was observed on the liner 
surface even after repeated runs. A check of exit temperatures at five locations 
in the e^diaust for (F/A)q^ « 0.016 (air flow > 14.2 M^/Min) showed the center- 
line temperature to be about 500°K above that for the original liner. Exhaust 
temperatures dropped quickly, however, with distance from the centerline and 
the resulting mean temperature was not significantly higher than that for the 
original liner. 
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Abstract 


Experinents and analysis on a gas turbine coabus* ** *** 
tor unit are reported with a view in Bind to sepa- 
rate propagated acoustic power from non-propagating 
"pseudo-sound". Analytically, it is suggested that 
a transition frequency will exist below tdiich the 
interior pressure fluctuati<ms are non-propagating, 
whereas above this frequency, of the order of 100 Hx, 
the noise is doBinated by propagating acoustic waves. 
Coherence aeasureaents are reported tdiich show this 
corc..pt to be borne out experiaentally. Coherence 
between interior and exterior nicrophones is sBas- 
ured over a wide range of experiaental conditions 
for a gas turbine coriwistor esdtausting directly to 
the atBosphere. The ;wrposc is to show that below a 
certain frequency, measurenents of interior noise are 
not indicative of coBbustion noise •iltimately prop- 
agating froB an engine. 

Honenclature 


cosibustor radius 
speed of sound 

denoBinator defined by Eq. (15) 
specific sensible internal energy 
transverse tail of series defined by Eq. (17) 
defined by Eq. (11) 

cross power spectrum between signals i and J. 

, ,. 1/2 


Bessel function of m order 

wavenuad>er defined by Eqs. (13) 
length of combustor 
Hach number 


Q heat release rate per unit volume 

r radial coordinate 

5 cross-section area 

t time 

u axial velocity cenponent 

V can volusB 

v vector velocity 

X axial coordinate sKsasured from the upstream 

end of the combustor 
3 specific acoustic admittance 

X eigenvalue defined by Eqs. (IS) 

C constant defined by Eqs. (IS) 

Y ratio of specific heats 

A constant defined by Eqs. (15) 

¥ transverse eigenfunction defined by Eqs. (15) 

T) defined by Eqs, (IS) 

p density 

6 angular coordinate 

111 frequency 

] of the order of 
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Subscripts 

c acoustic (dilatational) quantity 

e casdMistor exit plane 

mn suaaution indices for transverse modes 

o dumay integration variable 

v vortical (turbulent) quantity 

w side walls of Che combustor 

ai Fourier transform 

Superscripts 

~ dimensional quantity 

- time sman quantities or ensemble average of 

statistical quantities 
a sunaCion index for transverse modes 

' perturbation quantities 

Introduction 

It is reasonably well understood that in order 
for noise generated interior Co a curbopropulsion 
system to propagate to the surroundings a dilatation 
wave must accompany this noise. That is, in order for 
a propagacional wave to be present, which carries 
acvjuscic energy, the divergence of the velocity vec- 
tor Bust be non-zero. Alternatively, saBe of the in- 
terior noise is caused by vortical Butions cf the 
turbulence which are local in nature, not propaga- 
tional. This problem of "pseudo-sound" coexisting 
with propagacional noise causes difficulties in in- 
terpreting interior microphone laeasureniencs insofar 
as their connection with Che ultimately radiated 
sound is concerned. Nevertheless, one of the priaiary 
techniques used to Investigate the existence of 
"core engine noise" or "excess noise" has been the 
comparison of interior and exterior microphone meas- 
urements (1"^) 

A technique used previously^' to detect a 
causal relationship between an Interior and exterior 
microphone has been the cross-correlacitm function. 
The use of this function in interpretation of the 
results requires exact specification of Che two 
functions which are directly proportional to each 
other. So, for example, the interior pressure and 
the tiBB derivative of the pressure were cross- 
correlated with the far field pressure in Ref. (3) 
to determine which of the two functions was related 
in direct proportion to the far field pressure. A 
less specific but powerful function to investigate 
a linear, causal relatlm between two variables is 
the coherence function.' ' It has the property that 
if two variables are related by a linear transform 
operation it will be unity, regardless of the trans- 
form operation. It will be zero if there is no 
causal relationship. This behavior occurs regard- 
less of the transform operation. For example, in 
Ref, 3 the coherence fund' between either the 
time derivative of the intei jr pressure or the pres- 
sure itself and thr fi’" tie pressure would be 
unity if either o/ le two ..antities were perfect- 
ly correlated with the far field pressure because 
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differentiation is a linear transform operation. The 
coherence function will be used in this paper to 
determine the frequency regimes, if any, in which 
there is a linear causal relation between interior 
and exterior pressure fluctuations. 

The purposes of this paper are two-fold. First, 
it will be analytically established that below a 
certain frequency (apparatus dependent) the hydro- 
dynamic noise generated by turbulence within a com- 
bustor is non-propagational and consequently inco- 
herent with exterior sound laeasurcments. Secondly, 
for a particular experimental set-up, this fact will 
be demonstrated. 

Analysis 

The theory of coaibustion noise for the gas tur- 
bine combustor will assume that the fluid mechanics 
process Is one of Che gas phase alone. Therefore, 
the complications of fuel evaporation and mass, mo- 
mentum and energy transfer between liquid and gas 
phases are avoided. Consideration of these processes 
is not necessary Cor an adequate understanding of 
combustion noise generaCloa, since the dominant 
noise generation process is a gas phase one - die 
unsteady heat release In the turbulent gss phase. 
Although wss release fluctuations in space and tism 
by evaporation frcn th« spray droplets would indeed 
cause some noise, or let of magnitude arguments sup- 
port Che contest ion that the extreme heat release 
involved causes greater local fluctuations in gas 
velocity, the primary source of noise. The approach 
will be CO study the behavior of Che fluctuating 
components of various field quantities about their 
Cime-oean values, these fluctuations will contain 
both Che turbulence and acoustic fluctuations, and 
a split into vortical turbulent motions and acoustic 
quantities is found possible by Che method below. 

Molecular transport processes are neglected In 
comparison with turbulent transport. Althou^ it 
would be necessary to consider these to recover ex- 
act flame structure (since the actual combustion 
process must take place in the microscale of Che 
turbulence), the flame structure is basically left 
as an unknown in Che following treatment, to be 
determined by knowledge of both Che time-mean field 
quantities and the fluctuation quantities. 

The equations of continuity, moraencum, energy 
and state for an inviscid perfect gas are 

^ ~ , — 'v ^ ^ ~ ~ 

rr" + V. (p V )“ 0, TC- + V . I? v » - -sf~ , 

Ot ■* ot 0 

+ "p ® ^ -fe- ■» ( y-1 )“e (1) 

D t 0 1 P P 

Ncndimenslonalizing Eqs. (1) by the can length, Che 
exit plane mean speed of sound, density, and pressure^ 
the conservation equations become 

+ 7. (C V ) . 0 (2) 


^ - V f- dT ■ ^ 

Constructing Che r.ime mean equation th usual re- 
sults are found. That is, the pressure iradienc is 
of the order of the square of Che Mach : isnber so 
that if M^_< <1, p fts 1 ; furthermore, Che primary 
effect of ^ is to cause a dilation of Che velocity 
vector, not a strong pressure change. These remarks 
presume well distributed combustion Chrou{^ouC Che 


coadnistor so that 7.v is of the order of unity 
throughout the combustor. This is generally true 
for a spiay type diffusion flame. 

How consider the linearized fluctuation equations. 
By linearization all processes In the turbulence re- 
lying on nonlinearities such as generation, spectral 
transfer and decay (which is not considered because 
viscosity is not considered) are omitted. Consldera- 
tion Is therefore focused on Che energy containing 
eddies over die can length which is presumed short 
compared with a length required for a substantial 
alteratlm of the cuthulcnce structures. Condilning 
Eqs. (2) and (4), Eqs. (3) and (4) became after 
linearization 

dv' 

^ + 7(v.w') + 7xvxv' + 7xv'x w • - 

at — - - — y p '5) 

1^ + v'.7p +V . Vp'+ vp'7. V + VP 7.v'» (v-l)Q' 

(6) 

In Eq. (6) the second term nay be dropped because 
7P Is^(H^) and in the fifth term p * 1 may be used, 
correct Co terms of^M^) compared to unity. Now let 
Che perturbation velocity vector be split into its 
lamellar and solenoids 1 compements. 


where 7, v ' « 7XV* » 0. Here v' is associated with 

V •'•SI 

a dilacacional acoustic motion and v' is associated 

-nr 

with vortical turbulent motions convecced by the 
mean fluid motion. Because of the linearity of Eqs. 
(5) and (6) Che pressure may be considered as split 
into two components p^ and p^ and Eqs. (3) and (6) 

may be split as follows; 


^ + 7 iv.vj ^7xvxv^ - 
dp' 

+ v,7p^ 


-r— + v.7p 
St -• 


— ^ j,' « - yp' 7.V 

Pt *^V .. 


r— - + 7(v,v )+7xvxv +7XV xv» it- 

3C P 

In Eq. (8) Q' has been associated with the acoustic 
problem since heat release per unit volirae can pro- 
duce no torque, only dilatation. In Eq. (9) the ex- 
pected result is present that the pressure fluctua- 
tions ptoduced by the vortical motion are convected 
by the man flow, modified only by the rl^t hand 
side of Eq. (9). If terms of the order of the Mach 
nmidier are neglected, compared with unity, the sec- 
ond and third terms of Eqs. (7) and (8) may be neg- 
lected If ^ is^l) so chat wavelengths arc of the 
dc 

order of the combustor length. This assures that 
Che 7 operator is of order unity. In this case, 
caking the Fourier transform of Eqs. (7) and (8) 
and '■ombining 


+ 

“au 


at L 


(v-l) = 


Equation (11) is an inhomogeneous Helmholtz equation 
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driven by the fluctuating heat release. Furthermore, 
the source term Is the local Eulerlao time derlva- 
t lye of the heat release as has occurred previous- 
I'F, ^not dte total time derivative as in the formu- 
lation of Ref. (6). 


At lOH values of the frequency Eq. (11) is not 
accurate and the primary information comes from Eq. 
(8) as 

W.v'-(y-1)Q' (12) 


neglecting terms of 0(N) compared to unity. The 
pressure is then determined by Eq. (7) without the 
time derivative tern. Equations (11) and (12) are 
limiting forms of Eqs. (7) and (8> at w of ^1) and 
II) -*0, respectively. Integrating Eq. (12) over the 
cmabustor volume and assuming a high Impedance head 
end and walls of the coodMstor, the Fourier trans- 
form of the result yields 




(U) 


Eq. (13) states that the average normal velocity 
fluctuation at the can exit plane Is directly caused 
by the sum of the heat release fluctuations Interior 
to the can. The average velocity fluctuation at the 
can exit will radiate sound in accordance with the 
isipedance relatitm of Ref. (7).Inthe limit of low 
frequency and Just outside the can 


Ing system. This will have to be determined experi- 
mentally. The transition will be called the hydro- 
dynamic to acoustic transition since Eq. (11) has 
wave phenomena clearly associated with it while Eqs, 
(9) and (10) are cimtrolled by the turbulence field. 
Consider, then, the solution of Eq. (11). To solve 
Eq. (11) the acoustic behavior of the liner walls 
and exit plane must be known. From the work of Ref. 
(7) It Is known that the acoustic impedance of the 
exit plane should differ little from that of an un- 
flanged pipe, open to the Infinite surraundlngs.From 
the work of Ref, (8) the head end should behave 
basically as a hard wall. The side walls present a 
complex acoustic picture, however. It If known ex- 
perimentally that only weak can resonances appear, 
so that the side walls must be reasonably absorbent, 
acoustically. An acoustic admittance, 8 , will be 
assigned to the side walls which will b¥ left as a 
free parameter. The magnitude of 8 can be deduced 
approximately by a comparison of ^e theoretical 
and experimental resonance peak heights in the pres- 
sure. Under the approximattixis leading to Eq. (11), 
therefore, the boundary csmditlons on p are (assum- 
ing a cyllncrica. cosibustor) 


X « 0 


3p 

-Is 

dx 


r “ a 


dp 

*ui 

dr 


+ i (I) 8 P “0 

W ID 


potiofv dS«®l^dV 
tu <1 <1 ® 

5e V 

p -• 0 as » -• 0 

09 

Consequently, since must also go :o zero If ij — 0 

(by definition of a fluctuating quantity) an observed 
experimental fact of very low exterior pressure fluc- 
tuations at low frequency is recovered analytically. 
However, it will be seen experimentally that at low 
frequency Che Interior pressure fluctuations remain 
high while the exit plane fluctuations disappear. 

An interior microphone will see both p^ and 
whereas exterior mlcroptores will <raly see the ef- 
fects of p,', it being a propagatlonal wave rontain- 
ing fluid dilatation. O') can therefore expect a 
possibility that below a certain frequency, if p^' 
dominates p^' at hl^er frequency, p^' may dominite 
Pg'. This is only a possibility analytically because 
both pi^' and p^' must go to zero as ®->0 but the end 
plane impedance relation might drive Eo zero 
faster than Any dominance of at any fre- 
quency naist come from the source Q'and at ® of0^1) 
Eq. (11) is the dominant equation If the source is 
strong enou^. Experimentally It is known that prop- 
agational sound, called combustion noise, is heard 
from combustors, and, consequently, there Is a fre- 
quency regime where Eq. (11) Is valid. The point is 
that the limit of low frequency there is no rea - 
son to expect a correlation between an interior 
pressure masiirCTent and t^ pressure measurement 
o utside the can . At sufficiently low frequency there 
is the possibility that non-propagating "pseudo- 
sound" will dqmlnate Interior pressure measurements. 

Because the exact nature of the turbulence field 
Is not known, it Is inposslble tu state the exact 
value of ® at ^Ich a transition will occur frnii 
Eqs, (9) and (10) to Eqs. (7) and (8 ) as the govem- 


3p 

0) 

3x 


+ i u; 8 P 
"e "id 


(14) 


«diere a is the can radius and 8 * u o c 

e ii)g e e 

the specific acoustic admittance of the exit plane. 


The solution to Eq. (11) subject to Eqs. (lA) 
may be worked out by standard methods, by expansion 
in terms of normal modes of the combustor. Since Q 
is a randcm function of position in the combustor 
can there are no sycmetries that may be invoked to 
simplify the solution and the solution will contain 
both sysmetric and anti-synmetrlc modes. Since in- 
terest is ultimately in oomparing a wall-measured 
pressure against the exterior radiated pressure, 
the solution Is here given for a wall pressure. The 
solution Is 


(x,a,e) 


m, n • 0 

a • + 1, -1 


+ P 

U , (1) , 

mn,l mn,Z 
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a "m L e t tan J 

‘’m “ i S k A D 

tmipX snn mn nin 


X a 

f dx cos k x ; dS f (x ,r ,0 )Y_(r ,6) 
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m + (B„ IS a ) 
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«> 

mn 


J-^m 


e = lro = 0,e = l/ 2 m >0 (15) 

Ihe **101 solution contains the contri- 

bution of f from X ° 0 to X and p contains the 

IB ‘^U )2 

contributions from x (the measuring point) to l.lhe 
first tern in the overall sum (m«*n =0 , o “ + 1 ) is 
the (nearly) plane wave mode. It is not absolutely 
plane because Koo is calculated for small 

[ Kqq a I from the Bessel function equation to be 


2 ia,S , 1/2 
K = ( !!_) 1 IK 

0 0 N / L o 




For practical purposes, however, little error is 
made in the use of Eq* (16) for IK ! as large as 

ooa* 

unity. The finite magnitude of K means that 

oo 


Y = J (K r) 1 

OO o oo 


X and 1, p - 0. Additionally, if f is nonsero 

IB , U) 

OO 

for only a short axial distance near x ■ 0 and m is 
low, a careful Inspection of Eqs. (IS) shows the 
anticipated behavior that all transverse modes ate 
cutoff (Fi^) if the measuring point is sufficiently 
far downstream. This latter case cannot be expected 
to occur in practice since at best one half of the 
can length must contain active combustion. On the 
other hand, there is a reasonable expectation that 
near the quarter wave resonance point the plane 
wave mode, p , will dominate F. Further- 

“ 00,1 and 2 

more, in the limit of small u) (m^O) the plane wave 
mode will dominate F if is small because of the 
appearance of (b in the denominator. Oo the 
other hand, it should be recalled that at sufficient- 
ly low IS Eq. (11) ceases to be valid. Consequently, 
numerical calculations are in order to determine 
Che actual magnitudes of various terms. 

In order to experimentally determine whether or 
not any terms may be neglected in Eo. (17), howevei; 
the coherence function may be used.l^) Consider a 
measurement of the pressure at two combustor itosl- 
Cion, X and 1. The cross power spectnmi of these 
two pressures is 

G s p (x) p* ( 1 ) 

Xl IS IS 

where the * denotes a complex conjugate. Here 

P|s ' Pis^ (1) +F (1) 

00 


P|B^*> “ P|S^ ^ pi 

00 ,1 00 

Now note that pact of p ^ (1) contains p (x)in the 
'** 00,1 '** 00,1 

integration of Eqs. ( 15 ). In fact p ^ (l)=H(s)p ' (x) 

IS . JJ , 

00,1 00,1 

-I- J where H(is) Is a linear determined transfer func- 
tion and J is a frequency dependent randem function 
which is incoherent with any other term. Consequently, 
let p ^ (x) s b and the pressures take the form 


has a radial dependence, and, consequently, p has 

^00 

a radial dependence which keeps it from purely plane 
behavior. 

For bg large, as it is at low is, -1, and 

Dqq has a minimum near is = n/ 2 , the quarter wave 
resonance point. Thus, a maximum occurs in 8 t 

the quarter wave frequency, as expected. Because 
each term in the stun contains a different functional 
weighting of f^ in the integrations, each term in 
the sum is partially incoherent in phase with the 
other, unless, of course, f^^ is not random in posi- 
tion, which is not to be expected, p^^ may also be 
represented by 


p^(l) » Hb + K 
Pj^(x) * b + L 


where K and L are incoherent with b. If a many sam- 
ple average is i aken of G , to yield G there re- 
sults ‘ ** 

G , = b b* H* 


The auto spectra after averaging are 


b b* + L L* 


p “ p + p + F 

IB IS , IB 9 

00 00 


where all three quantities are incoherent with each 
ocher and F is Che entire transverse mode call of 
the doubly infinite series. It would be nice if it 
could be shown that some of the terms of Eq. (17) 
could be neglected. Clearly, if f^. •> 0 between 


Gjj - b b* H H* + O*” 

The coherence function is defined by 


G G*, 
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and diete results 

2 _ H H* 

(bT* + n> )( rb* H H* + Flc* ) (18) 

If etcher a) events between x and 1 doninace be* 
havlor (J and L are large so chat K and L are large 
coapared with b , or b) F is dominant and F(l) Is 
incoherent with F(x) , is small. On the other 
hand if b is the dominant term ~* 1. A coherence 
function of unity (or near unity) in a given fre- 
quency range will therefore show that a single term 
dominates the solution of fiqs. (IS) and that events 
between x and 1 are unimoortant (^nO). Further- 
more, the cohereoec fimctlon will also show where 
Eqs. (7) and (8) take over from Eqs. (9) and (10). 
This can be seen if the microphone at x is dominated 
by and the x > 1 microphone is placed slightly 
outside of the cagmsCor end plane to sense only 
propagational sound, p'. 

If a frequency band is found where y n 1 • It is 
asstmed that f 0 between x and 1. Any deviation 
from y^ M 1 at'*^igher frequencies oust then indicate 
the transverse Diodes competing with the plane wave 
mode. Calculations have been performed on Eqs. (IS) 
to compare the magnitudes of p ^(x) and p'^^(x). A g 

tti ^10 ^ 

magnitude has been chosen to gi^ roughly »e quar- 
ter wave resonaice peak aiagnitude which is experi- 
laentally observed, has been assumed real, on the 
basis of no other information. In the calculation, 

T ^ (a, 6) ° 1 and K^g real (gwj'^3) have been assumed 

for the order of magnitude calculations. Ihe posi- 
ti<m X has been chosen, to correspond with the ex- 
perimental setup. Ihe expression, for is the low ... 
frequency asymptotic expression of Levin^Schwinger'^' 


ge “ l/[oo a( 0.6 i + 0.2S is a) ] 

This has also been used for the first transverse 
mode calculation. The elimination of the wall damp- 
ing in the transverse mode calculation will raise its 
magnitude compared with the case ot damping; conse- 
quently, this is a conservative comparison of the 
dominance of the plane wave mode over the transverse 
mode. Ihe integrals in Eqs. (IS) have been estimated 
by a) carrying out the 3 ^ integration as though 
were Independent of md b) the transverse inte- 
grations have been estimated as 



where A^ is an average (unknown) source strength. 
Shown in Fig. 1 are calculations for the plane wave 
mode and the first antisynmetric transverse mode 
(1,0), which is the first transverse mode to cut on 
as frequency rises. Starting from top to bottem, the 
plane wave mode at x ° 0.6S3 with source contribu- 
tions from 0 to X (00,1) is the dominant term. It 
shows the quarter wave resonance near uj ° tt/2. 
Sli^tly lower in aiagnitude is the plane wave mode 
evaluated at x ■ 1, now using the source contribution 
from X • 0 to x ■ 1. The reason it is slightly lower 
in magnitude, even though the source contributes 
over the full can length, is that Che Impedance con- 
dition at the can exit plane requires nearly con- 
stant pressure, whereas there is no such restriction 
interior to the combustor. Much lower in magnitude 
is the plane wave contribution to Che x 0.655 


poeltlou from the combustor position downststream of 
the transducer positlm. (00,2). Consequently, even 
if Che combustor source strength were as large down- 
stream as upstream of the transducer position it 
would not contribute to Che transducer oieasuremsnt. 
Finally, well down from the plane wave modes, is Che 
first transverse mode conCributton. It is almost 
completely cut off at the transducer location over 
Che frequency range shown. Theoretically, It would 
cut on completely at a dimensionless frequency of 
8.8 The frequencies of interest experimentally are 
well below this value so that the F function of Eq. 
(17) may indeed be neglected. The major conclusion, 
therefore, is diet in Eq. (18) ~ 1 should result 

over the frequency range of interest as long as the 
frequency is high enou^ that Eqs. (IS) are a valid 
solution to the problem. 


Experiment 


Apoaratue 

The gas turbine coodnistor used in this study was 
described in Ref. (10). A schematic of the setup is 
shown in Fig. 2. A water-cooled pressure transducer 
capable of extracting signal above 135 dB re 2xl0~^ 
N/m^ was located on the liner wall, 92 tag from the 
can exit plane. A near field 1/2 inch condensor 
microphone was located in the can exit plane, 113 mn 
to the side of the liner lip. This transducer meas- 
ures the can exit plane dynamic pressure as long as 
the frequencies are resr.ricted to wavelengths long 
compared with this separation distance. Far field 
microphones are employed to detect the onset of 
scattering from the burner hardware (noo'^onopole 
bdiavior of the can exit plane) as will be explained 
later. 

All signals are rec->rded on FM Cape for later 
digital Fourier analysis. The primary items of in- 
terest are spectral shapes aud coherence analysis. 
Calibration procedures are therefore unimportant, 
except to be satisfied that the microphones are be- 
having linearly in voltage output vs dynaiic pres- 
sure level. 

Results 

Figure 3 displays tne spectra obtained by fixed 
bandwidth frequency analysis of the interior micro- 
phone signal. Shown are Che runs at extreme ends of 
Che airflow and fuel/air ratio matrix. In all cases 
there are two regimes to the spectra. Below about 
200 Hz Che spectra arc quite flat. Above this value 
Che spectra undergo a hump. The upper frequency 
limit was chosen as 800 Hz, for reasons to be appar- 
ent later. Depending upon fuel/air ratio (speed of 
sound) a weak quarter wave resonance is seen in all 
curves in the vicinity of 400 Hz. 

In contrast, look at Figure 4 which displays a 
typical near field microphone spectrum (although 
here a wider frequency range has been considered). 

Ac low frequency the speernim falls off at roughly 
13 dB/octave whereas the interior microphone oieas- 
ured a flet spectrum. By way of interest Figure 4 
shows the 1/4, 3/4 and 5/4 wave resonance peaks. 

The major point is, however, chat the near field 
spectrum has behavior similar to the interior mlco- 
I^one spectrum only above about 200 Hz. More docu- 
mepCation on this effect is shown by the several 
near field spectra in Fig. 5. Also shmn is a typi- 
cal far field microphone spectrum which is similar 
in shape to Che near field spectra, except that 
ground reflection effects are seen in the spectrum. 
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Now consider Fig. 6 «Alch shows Idle near to far 
field coberenc'* which is hl^ (nearly perfect co- 
herence) over the range 100-600 He and Is acceptably 
hlgb even below 100 He. As long as the frequency Is 
suff i.clently law tdiat the Is behaving ns a sumo- 
pole radiator the far field Is quite coherent with 
the near field even In the presence of ground ref lec- 
ions; this occurs because ground reflections cause 
only a linear transform operation on the near field 
signal to recover the far field signal. Furthermore, 
In the freipiency domain the distance traveled to the 
far field microphone only Introduces a phase lag, 
again a linear transform operation which cancels out 
In coherence analysis. Above 600 Us one concludes 
that scattering from the burner hardware is such 
that Bumopole behavior Is destroyed; a single point 
sieasurement In the near field is Insufficient to de- 
scribe the total characteristics of the far field 
radiation. Con>eq>'eatly, ,-elow 600 He a near field 
measurement Is adequate to describe the spectrim of 
the ultimately radiated soud power, and attention Is 
focused only below 600 Hz. 

Now view Fig. 7 which portrays the coherence be- 
tween the Interior and near field microphones for 
several run conditions. Hie coherence behavior Is 
Che same for all runs; It Is above O.S between about 
150 He and 600 Hz and Is nearly unity over the range 
200-500 He. As frequency drops below 150 Hz the co- 
herence drops to nearly ;;uro. Hie interpretation of 
these results Is as follows, a) a source other than 
coeoustlon noise Is daolnacing the Interior micro- 
phone measurement at frequencies below ISO Hz, b) 
be:muse of the above analysis and the non-corres'- 
pondence of the spectral shapes of the two micro- 
phones at low frequency the othet noise source Is 
hydrodynamic noise caused by the vortical curlmlent 
motions and those motions do not provide propaga- 
tional sound and c) above 600 Hz Che low coherence 
Is caused by the failure of a single point measure- 
rent In the near field Co characterize the entire 
radiation pattern; tiie monopole radiation pattern 
has ceased. Hie drop In the near field spectrum at 
low frequency Is faster than a dropoff proporclcnal 
to frequency alone, indicating that is falling 
off at roughly 6-7 dB/octave below about 300 Hz. 

Hils appears to be Che primary reason for the take- 
over of hydrodynamic pseudo-sound at low frequency. 

Discussion 

Hie results presented here can clearly be expect- 
ed Co be apparatus-dependent, especially with regard 
to the frequency at which hydrodynamic noise and 
combustion noise are equivalent at an interior loca- 
tion, The only contention here is that generally 
such a transition at some frequency may be expected. 
In Che experiment here there were only two basic 
noise sources - flow noise and combustion noise. In 
an Installed configuration at least two other com- 
bustor-turbine interaction noises will come into Che 
picture to confuse the results; there will be entropy 
nolse^^ and vorticity-nozzle interaction noise, 
both of which will cause propagational noise and 
which will be sensed by both interior and exterior 
microphones. These sources have been avoided here 
by directly discharging to the atmosphere. 


rations Is required to dtceiriue the ariount of 
pseudo-sound In the ccmbuscor as compared with prop- 
agaClonal sounc. 

In die experimeuCs here It was fortunate chat 
the frequency content of combustion noise was low 
enough and the burner else was small enou^ chat 
a) Che plane wave mode dominated the Interlpr acous- 
tics and b) the can radiated basically as a monopole 
over the frequency range Chat cmtalned most of the 
combustion noise. Severe complications in Interpre- 
tation can be expected if this is not the cas^. 

Reference ;: 

1. Abdelhamld, A. N. , Harrje, D. X., Piece, E. G. 
and Sunmerfleld, M., “Noise Characteristics of 
Cadnisclon Augmented Jets at Midsubsonlc Speeds," 
AIAA J. 12, 336-392 (1974). 

2. Gerend, B, P., Kumasaka, H. A. and Boundhlll, J. 
P., "Core Engine Noise," Aeroacoustlcs : Jet and 
Combustion Noise ; Duct Acoustics (Nagamatsu,ed.) 
AIAA, New York (1975) pp. 305-326. 

3. Fleet, E, C., Teshner, M. D. and Susmerfield, H. , 
"Combustion CeooMCry and Ccmbustltm Roughness 
Relation to Noise Generation," Second Interagon- 
cy Symposium on University Research in Transpor- 
tation Noise, Vol, 11 (1974) pp. 723-729, 

4. Bendat, J. S., and Plersol A. w. , Random Data : 
Analysis and Measurement Techniques . le/- 
Incersclence, New York (1971) p. 32. 

5. Strdtle, W. C. and Shivashankara, B. N. , "A 
Rational Correlation of Combustion Noise Results 
from Open Turbulent Flames," Fifteenth Symposium 
( Incematianal ) on CosJiustion . The Combustion 
Institute, Pittsburgh (1975) pp. 1379-1388. 

6. Chiu, H. H. and Summer fie Id, M., "Theory of 
Combustiun Noise," Acta Astronautics . 1, 967-984 
(1974). 

7. Hahbah, M. M. and Strahle, H. C., "A New Method 
for Solving Problems of Sound Radiation for t'ue 
Case of Low Mach Number," AIAA Paper No. 76-582 
(1976). 

8. Janardan, B. A., Daniel, B. R. and Zinn, B. T. , 
"Characteristics of Response Factors of Coaxial 
Gaseous Rocket Injectors," NASA CR-134 788 0975). 

9. Levine, H. and Schwinger, J., "On the Radiation 
of Sound from an Unflanged Circular Pipe," 
Physical Review . 37, 383-406 (1948). 

10. Strahle, H. C. and Shivashankara, R. N. , "Com- 
bustion Generated Noise in Gas Turbine Combus- 
tors," Engineering for Power . 98 , 242-246 ',1976) 

11. Cumpsty, N. A., "Excess Noise from Gas Turbine 
Exhausts," ASME Paper 75-6T-61 (1975). 

12. Pickett, G. F. , "Core Engine Noise due to Tem- 
perature Fluctuations Convecting througii Turbine 
Blade Rows,” Aeroacoustlcs : Jet Noise . Combustion 
and Core Engine Noise (Schwartz, ed. ) , AIAA, New 
York (1976) pp. 589-608. 


A major point to be stressed is that combustion 
noise may be present in the far field signature even 
if there is little coherence between internal and 
external microphones. This was precisely the case at 
low frequency in these experiments. Consequently, 
some care in experiment design in installed configu- 
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Figure 1. Relative magnitudes o£ various mode 
Strengths at various combustor posi- 
tions £or a given combustion noise 
source strength* 


B SK FAR FIELD 

MiCROPMONc ' 


I 


( 


SBK NEAR field 
KICRORHONC 

AVL PRESSURE ^ 
TRANSDUCER' 




AIR 

FU3« 



Figure 2. Schenactc of experimental setup. 
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Figure 3. Interior microphone spectra as func- 
tions of airflow rate and fuel/total 
air. The fraction of air bypassed is 
approximately 26%, 



Figure 4. Typical near field microphone 
spectrum. 
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Figure 7. Coherence between the Interior and 
exterior microphones for several 
run conditions. 


Figure 5. Near field spectra as a function of 
airflow and fuel/total air. Also 
shown Is a typical far field 
spectrum. 
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Appendix C 


Theory of Combustor Noise 


Formulation and Solution 

The theory of Appendix B Is used as the starting point here. The differen- 
tial equations are based upon the linearized vortlclty - acoustic field eqiia- 
tlons.^^^^ Equation (11) of Appendix B is the acoustic equation to be solved, 
«dilch Is here reproduced as 


P^ + « P. 




- f 


(tt 


(Cl) 


where the nomenclature is defined In Appendix B. Equation ( Cl ) Is an acoustic 
field equation, independent of the turbulent, vortical part of the flow (by 
assumption), except as the two fields may Interact at boundaries. Equation (Cl) 
was derived on the basis of a low Mach number flow. 

In Eq. (Cl) die acoustic pressure is driven by fluctuations In the heat 
release throughout the combustor, represented by f^. I3ie boundary conditions 
applied to Eq. (Cl) In Appendix B are given as Eqs. (14) In that appendix. 

The boundary condition at the combustor exit plane (x=l), however, was not the 
most general boundary condition, and here such a condition will be introduced. 
Since It was found In Appendix B that only the plane wave mode was Important 
at the experimental frequencies of interest, only the plane wave mode will be 
treated here. 

The most general boundary condition at the exit plane for the plane wave 
mode Involves the axial velocity, pressure and entropy fluctuations In the 
form<"2) 


u 




Pc 




+ or CT 


0 )- 


(C2) 


"e e 

where a Is the entropy admittance coefficient, and the dimensionless entropy 

f luct\iatlon Is defined by ^ ~ / 

o = -rsr 
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and Is equal to 


/ 

a 




i! 

p 


for perfect gases. The oscillations of Eq. (Cl) are isentropic (a'=0) and 
the primary source of o' is in temperature fluctuations produced by unequal 
heating of different fluid elranents as they travel down the combustor. Diis 
is emphasized by Eq. (4) of Appendix B which may be written 


Q 

Dt Y 

Both the acoustic pressure and the entropy fluctuation are consequently driven 

by heat release fluctuations. Moreover, they interact at the boundary through 

Eq, (C2), It must be emphasized, therefore, that entropy noise (noise due to a ) 

and direct combustion noise (noise due to f) are related, but not identical. 

The superscript on a in Eq, (C2) is to emphasize that the plane wave 

*“e 

mode requires a cross section average of o ; that is. 



(C3) 


Using the nondimensionalization scheme of Appendix B and the fact that a' is 
dominated by temperature fluctuations 




f 


where T is temperature. These temperature fluctuations are presumed known from 
some independent measurement. 

Ihe boundary conditions at the side wall and head end are the same as in 
Eqs. (14) of Appendix B. The solution to Eq. (Cl) then follows from standard 
methods for an interior wall pressure measurement point downstream of the 
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combustion zone 
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where the nomenclature is the same as in Appendix B. The effect of the exit 
plane entropy fluctuation Is to add a term to the solution proportional to 

A 

ao • In the case of experiments with no nozzle o = 0; hence there Is no 
“e 

effect of the entropy fluctuation In noise generation. 

Now consider three microphone positions at x = Xq, x»l and at the nozzle 

exit plane. Denote the plane wave pressures there as p , p , and p .In 

U)q’ *^u>2 

general the true measured pressure at an Interior microphone position will be 
contaminated with hydrodynamic noise which does not propagate to the near 
field, considered as the p^ microphone If it Is placed near, but outside of, 
the jet stream. At x»x^, Eq. (C4) may be rewritten as 


e 


(C5) 


where H Is the hydrodynamic noise, f^ and f^ are frequency dependent transfer 

functions and A is the transform of the weighted heat release fluctuation 
u) 

A = f cos k-rt X dV 
(U J U) 00 

V 

At the nozzle entrance plane, the acoustic pressure will also be computed 
from Eq. (C4) with an additional hydrodynamic noise contribution 



A^ f 2 ( UJ ) + ^4 ( “> ) + ( u) ) 


(C6) 


Finally, there are the acoustic transmission generation and reflection 
processes due to the presence of a nozzle at the exit plane. It will be 
assumed that the nozzle is sufficiently short (compared with a wavelength 
or length of an entropy pulse) that quasi-steady short nozzle theory may be 
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used.'' * In such a case the two conditions are a) conservation of mass and 

b) conservation of stagnation energy. For the acoustic part of p , which 

^a 

excludes these conditions yield 


Vpj I «2 / «1 
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“i s 




2 2 
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with z => 1 + ’.d^, and Isentroplc oscillations have been assumed. Since in 

the absence of entropy fluctxiatlons Eqs. (C7) may be expected to have a solu- 
tion, the determinant of the solution matrix must be zero, yielding 

(1 -Zj/Z2 ) + P 2 (1/^2 " h “2^^2^ 

^1 “ z (1 + 

^ — ^.Mi(l-f.p2/V 


Moreover, Eqs. (C7) must yield a solution with a finite. Consequently, 

'**e 




- (1 + 


Bj/Hj) Mj] 


Equation (C9) emphasizes the previously mentioned fact that as M 2 -*M^, a-*0. 
The only Item left to close the solution Is ^ 2 * Reference (C3) shows 
that at low values of M 2 the fact that there Is a temperature mismatch be- 
tween the combustor Jet flow and the surroundings does not greatly alter ^2 

(C4) 

over its Levlne-Schwlnger value.' ' Using the low frequency limiting form, 


’ 2 j ^ ^^2 0»25 ^^ 2 ] 


I 


I 


vhere is the radius of the nozzle outlet. On the other hand» near M2*>1» 
the choking condition yields ^2 * ^or the approximate calculations 
here a linear relation Is used for Mach numbers between 0 and 1; that is. 


p2 - ^2 (1 -M 2 ) + M 2 ( Y -D/2 

LS 


(Cll) 


Hie first of Eqs. (C7) and Eq, (C6) may then be written as 

r QfYao)g 1 . , 

\ “L\ Mj (1+P2/M2) J Y 

a 

“ j \ ^3 ‘^0,^ [ M^a+p2/M2> ^ 4 ] [ ^ 
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(C12) 


The expeilments of Appendix B verify that the hydrodynamic noise, and 
is not propagatlonal to position 2 and no hydrodynamic noise is present In 
the P 2 signal If the microphone Is slightly outside of the jet flow. Com- 
bustion noise dominates the situation, at least when no nozzle is present. 

Now the major fact revealed by Eqs. (C6) and (C12) is that entropy noise 
and combustion noise have different transmission characteristics to the P 2 
mio '>phone (which is also coherent with the far field microphone at the fre- 
quencies under consideration). From Eq, (C12), p^ depends upon both p^^ and 
a , but from Eq. (C6) p depends on A (combustion noise) and a (entropy 


U)i 


U>1 


U) 




noise). Depending upon the behavior of f^, f^, f^, and f^, which are all dif- 
ferent, It Is possible to have the following possibilities: 

Position 1, a dominant. Position 2, a dominant; Position l,o dominant. Position 
2, p dominant which means a Is still dominant; Position 1, A dominant. Position 
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2»a dominant, which means that combustion noise dominates the combustor 
Interior while entropy noise Is the strongest exterior noise; Position 1, 

A dominant. Position 2, p dominant which means that combustion noise 
dominates the exterior noise* 

Hie distinction between an entropy noise and a combustion noise domi- 
nated case becomes a bit unclear when one recalls that they are both funda- 
mentally linked to the heat release fluctuations. For the moment, however, 
consider that A and a are statistically uncorrelated, and this point 

U) UJg 

will be returned to In a later appendix. 


Calculations 

Assuming statistical independence of A^, and H, the ensemble average, 
one sided spectral densities of the measured pressures may be constructed from 
the Fourier transforms in Eqs. (C5, C6 and C12). Denoting a one sided spectral 
density by G with a subscript for the quantity under consideration, there re- 


sults 


°P„ ‘ “a * 's * Si, 


\ - «A b' <"3 b' 


G = G. U \ + Ge 

Po (Pi“H.) 5 s 6 


(C13) 


The spectral densities are dimensionless by construction when is the pressure 
used for nondimenslonsllzlng and the time unit is -t / c^. In the calculations 
all numbers have been plotted as a spectrum dimensionless in pressure units 
but on a per Hertz basis. 

Consider first the case with no nozzle. Sample calculations have been made 

3 

for the 19.8 m /min air flow case at an overall fuel/air ratio of 0,016; Figure 
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7 of Appendix B shows the coherence between the Interior and exterior micro- 
phones for this run. In the experiments the position of the interior micro- 
phone is at X ■ Xq ■ 0.655. For this case “ 0.1 and a « 0, so only 

and enter the picture. A spectral peak corresponding to the first quarter 
wave resonance occurs at 450 Hz for this case and experimentally there Is nearly 
perfect coherence between p^ and P 2 . Moreover, Independent experiments have 
shown that p^ senses no hydrodynamic noise at its location out of die jet flow. 
Consequently, at 450 Hz both p^ and P 2 are sensing combustion noise almost com- 
pletely. The experimental spectrum for p^ is shown in Fig. Cla. A value of 0 

U (I) 

2 

is chosen so that the theoretical spectrum (f^^ ) best osculates the spectral 
peak of Gp^; this value is 0.06 which indicates a hard, but not too hard, be- 
havior of the liner walls. This determines G^ at the 450 Hz point. Since t’.e 
physical location of the P 2 microphone is not precisely at the exit plane, a 
calibration factor to multiply the experimental Gp^ is obtained knowing G^ 
at 450 Hz and 

Sj - ^5^ 


I 1 

I • 


I 


Ihe calibration factor is then applied to the measured P 2 spectrum and the 

corrected spectrum represents the P 2 spectrum which would have been measured 

at the exit plane. This is shown in Fig. Cla. Since G senses only combustion 

2 

noise, G. may be calculated by G. = G /f_ as a function of frequency; this is 
A A J 

^ 2 
shown on Fig. Cla. Knowing G. , one can calculate G„ from G,, = G„ - G. f- . 

This is shown in Fig. Cla and is the derived hydrodynamic noise, lliis subtraction 

process becomes highly sensitive to small experimental errors when G^ falls more 

than 5-6 dB below G ; if this happens GL is arbitrarily set to zero. One sees 

that the interior microphone is dominated by hydrodynamic noise at frequencies 

less than, say 150 Hz. 


. 0 *: * 
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As a check on this procedure, ^Ich assumes only two noise sources present 

at the interior microidione and one source at die exterior microphone, the co- 

(C51 

herence function between the p^ and microphones may be computed. This 
is 


1 + 


^0 


It is comiHited and compared with experiment in Fig. Clb and the agreement Is 
very good. Ihe computed function is hi^er than the experimental function 
everywhere, which it should be because in reality there is some electronic 
noise and a low level jet noise which contaminates the two signals. 

Guided by this successful separation of coinbustion and hydrodynamic noise, 
it now remains to calculate the effects of entropy noise. For this, independ- 
ent measurements are required for . These ar: ^escribed in Appendix D. For 
these calculations Gg is known from those experiments. It is assumed in the 
calculation that remains the same with the nozzle on or off. Calculations 
are made for two cases, one with a contraction ratio of 3:1, which at the same 
mass flow and fuel/air ratio yields = 0.1, = 0.33, and one for a contrac- 

tion ratio yielding M 2 = 1.0. It is moreover assumed that is unchanged with 
a nozzle on (since the pressure is only varied over a factor of two in the com- 
bustor, this Is not expected to yiela a large error). 

Shown in Fig. C2a are the components of the calculated G , where GL has 

Po ^0 

been omitted for clarity (see Fig. Cla), for the "short nozzle" case (low con- 
traction ratio). It is seen that the entropy noise contribution to the signal 
is small, so that for practical purposes the interior microphone would be dom- 
inated by combustion noise. Notice that there is a strong change in the com- 
puted spectral shape, as compared with G in Fig. Cla. Since the calculations 

Pq 

were only carried to 600 Hz the resonance peak, which has shifted to a higher 
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field microphones for the no-nozzle case 
theory. 
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frequency in the presence of a nozzle , does not appears Continuing, G 


t 


(Pi-Hjl) 


is computed and then G . The results are shown in Fig. C2b. It is evident 

P2 

that the entropy contribution is still negligible, but approaching the com* 
bustion noise curve at the higher frequencies. 

The calculations are reported for the choked nozzle case and are shown 
in Figs. C3a and C3b. The results are that the interior microphone is domi* 
nated by combustion noise but the transsJ.tted noise becomes entropy dominated 
above 210 Hz. Moreover, there is an elevation of the high frequency content 
of the noise. 

Conclusion 

The conclusion is therefore clear, but may be due to the particular hard- 
ware used in the experiments. At sufficiently hi^ frequency and under choked 
or nearly choked conditions, entropy noise may be the dominant noise source 
in combustors. Under sufficiently low nozzle pressure drop combustion noise 
will dominate. Hiese theoretical results are borne out by the experimental 
results in Appendix D. It must be cautioned that because of the many approxi- 
mations used in the analysis that the results should be viewed in a qualita- 
tive sense only. Moreover, the assumption in the calculations that G. and G„ 

A 0 

are independent is not strictly true, as will be seen in Appendix E. 






Figure C3a. Choked nozzle interior mlcro(^one spectra 
and the combustion and entropy noise com- 
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Appendix D 
Entropy Noise 

Genera l 

This section deals with combustion and entropy noise radiation characteristics 
and the experimental procedures for separating them. Combustion noise radiation 
characteristics and the ways of separating them from the hydrodynamic noise 
have been explained In detail In Appendix B. Entropy noise Is generated In tur- 
bopropulslon systems when combustion generated hot spots are accelerated throu^ 
a region of mean pressure gradient. In the present program, such a pressure 
gradient Is Imposed on the hot gases by attaching to the burner exit a conver> 
gent nozzle of area reduction 3:1 or an orflce plate with 9 holes of 14 mm di- 
ameter each and having an effective area reduction of 10.5:1. An exit Mach num- 
ber, M23 0.20, Is obtained with the above nozzle whereas at hl^ Mach numbers, 

M 2 = 0.6, 0.8 and 1.0, the jet Is exhausted through a multi-hole orflce plate, 
rather than a single opening nozzle, to keep the jet noise Interference low. 

Experimental Apparatus 

The facility used for the present investigation Is the same as that described 
in Appendix B. The experlmentr.1 set-up with the relative locations of the trans- 
ducers is shown In Fig. Dl. The pressure fluctuations inside the combustor and 
the near field and the far field radiated sound pressures have been measured 
through a photocon pressure transducer and Bi^el and Kjaer type 4134, 12 mm 
condenser microphones. The reason for Including a near field microphone in the 
measurement scheme is twofold. At higher exit Mach numbers, the jet noise is 
expected to contaminate the exterior radiated core noise signals. In that case, 
the near field microphone by virtue of Its location, will receive more contri- 
bution from core noise radiation with a lesser degree of contamination from a 
small local area of the jet. However, the far field microphone signal may be 
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Figure Dl. ‘?chematic of experimental set-up. 
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donlnated by Jet noise as It sees the vhole Jet. This means that the near field 
signal, rather than the far field one, can be used with more confidence to rep- 
resent the exterior radiated core engine noise. Die second reason is that the far 
field microphone spectrum is masked by severe ground reflection effects, as shovn 
earlier in Fig. 5 of Appendix B, \diereas the near field microj^one, by virtue of 
its close location to the burner exit, is free from these effects. Since the en- 
tropy noise study requires the details of the combustion generated hot spots, a 
measurement program has been carried out to obtain the temperature fluctuation 
characteristics. Fast response chrome 1 -a lumel thermocouples of wire diameter 0.025 
tom are used to record the temperature fluctuations at the burner exit plane as 
shown in Fig. Dl. The time constants of the thermocouples are determined as ex- 
plained in Appendix F, and the tliermocouple signals are compensated for these 
time constants later during data analysis by suitable modification of the Fourier 
Analyser programming. All the signals, after amplification, arc recorded on an 
Ampex FR 1300, 14 chani.el magnetic tape recorder at a tape speed of 30 Ips. The 
recorded signals are then subjected to spectral analysis and coherence function 
estimates through a HP 5451A Fourier Analyser system. 

Results 

Spectral Characteristics : An exit Mach number, M2“0.05, represents the burner 
can operating without a nozzle while M 2 » 0.20 belongs to the convergent nozzle case 
and M2**0.6, 0.8 and 1.0 come through the use of the orifice plate. Hie airflow rates 

3 

vary between 7.0 and 19.8 m /min., with the overall fuel/air ratio maintained constant 
at 0.016 for all cases. Out of a wide range of test runs, only the above conditions 
are selected for the purpose of comparison of the results. Fig. D2 shows the spectra 
of the Interior pressure fluctuations. The upper frequency limit is chosen as 2000Hz 
because the information of interest falls below this limit. Figure D2a illustrates 
the spectral changes with the increasing pressure drop across the combustor as the 
exit termination is changed for a fixed airflow rate and fuel/air ratio, whereas 
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Figure D2, Interior photocon spectra as a function of 
airflow rate and fuel /total air. 



Figure D2b shows the spectral behavior with the increasing pressure drop across 
the orifice plate for varying airflow rates and a fixed fuel/air ratio. It can 
be seen from Fig. D2a that as the exit contraction ratio is increased, the spec- 
tral level rises significantly with an Increase in the exit Mach number, in the 
lower frequency range of 0-300 Hs. Around 300 Hz, all attain almost the same 
level, and then they fall off rapidly at tite same rate (10 dB/octave) with in- 
creasing frequency. Another observation made here is that the Interior spectrum 
corresponding to the exit Mach mmiber, M^^O.OS, eidiiblts a quarter wave reson- 
ance peak around 500 Hz. With an Increase in the exit contraction ratio corres- 
ponding to hl^ Mach number cases, this resonance peak should shift to a higher 
frequency range (800-1000 Hz), tending towards a half wave resonance. On the 
contrary, these peaks have disappeared in the high Mach number cases. Such peaks 
have been seen in experiments not reported here, when using an AVL interior trans- 
ducer. The occurance of a peak around 100 Hz in the case of high Mach numbers is 
explained at the end of this section. 


Fig. D3 shows the near field spectra for the same run conditions. While at 
low exit Mach numbers, there is a familiar combustion noise hump, high Mach num- 
ber cases eidiibit no such characteristics. Ihe near field spectra, pertaining to 
hi^ Mach number cases, do not resemble the usual combustion noise spectra. Hie 
ch£Uige in the spectral shapes at high Mach numbers could mean either entropy 
noise domination or Jet noise contamination of the near field signal. To investi- 
gate the above speculation, two checks were made. First, the scaling of the near 
field sound power level with exit velocity was obtained. The exponent on velocity 
in high Mach mmiber cases was higher than value obtained at low Mach number cases. 
However, the value of this exponent was far below the one corresponding to the 
jet noise velocity scaling law. This means that the near field microphone receives 
a greater contribution from core engine noise radiation although there is jet 
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noise con::amlnatlon to a lesser degree. Secondly, a spectral analysis of the 
near field signals was carried out. Fig. D4 shows a comparison of the near 
field spectra of the hot eidiaust at M 2 » 0.6 with the cold flow of identical 
exit velocity. It can be seen that tlie spectrum level with the combustion on 
is higher than that of the exit-velocity -matched cold flow by at least 3 to 4 db 
in the frequency range of 150-1000 Hz. This gives some confidence that the core 
engine noise dominates the jet noise, at least in the frequency range of interest 
to the present investigation. Later in this section, it is shown through coherence 
function analysis that it is indeed the entropy noise that contributes signifi- 
cantly to the near ^ield signal at high exit Mach numbers. 

Typical spectra of the temperature fluctuations the burner exit plane are 
shown in Fig. D5 The spectral shapes of this temperature fluctuations, corres- 
ponding to different exit terminations, as seen from Fig. D-5, are found to be 
very similar. Ihe temperature fluctuation spectrum is a broad band one and low 
frequency in nature. The spectrum gradually falls off with an increase in fre- 
quency. The r.m.s. temperature fluctuations are found to vary within a range of 
9 to 14% of the mean temperature for various cases. 

The theoretical evaluation of entropy noise requires a knowledge of the 
spatial correlation length scales of the temperature fluctuations. Two thermo- 
couples of the same size (0.025 mm wire diameter) have been used for this correla- 
tion study. One of the thermocouples monitors the temperature fluctuations at a 
fixed place at the burner exit plane while the other one registers the signals 
at various locations, moving radially away from the former. The fixed and the 
moving thermocouple signals have been cross-correlated and normalized to obtain 
the cross-correlation coefficients. The cross-correlation coefficient, as a 
function ot the separation distance between the two thermocouples, is shown in 
Fig. D6 for two cases corresponding to the thermocouple traverse along two 
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Figure D5. Temperature fluctuation spectra for 
typical operating conditions. 
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different radial directions. In both cases, the temperature fluctuations seem 
to be correlated over a distance of about 12 mm at the burner exit plane, thereby 
suggesting a temperature eddy of about 12 mm In size. The constancy of the corre- 
lation length scales In various radial directions demonstrates the axl-symmetric 
nature of the thermal eddies. Hie similarity in the correlation analysis results 
for two different test conditions is brought out In Fig. D? In all the above 
cases, the cross-correl^itlon is found to be positive everyxdiere, with a long tall. 

It is important to note that the spectral characteristics of the temperature 
fluctuations mentioned in the above paragraphs were obtained through single point 
thermocouple measurements at the burner exit plane. However, the entropy noise 
calculations require the quantities averaged over the burner exit area. The re- 
lationship between the single point thermocouple measurements and the area averaged 
quantities arc established below. 

Let t' and T represent the fluctuating and mean temperature measured at a 
point at the burner exit plane by a thermocouple and a ratio be defined as 




(D-1) 


where s' is the fluctuating entropy component and Is specific heat at constant 


pressure. However, the entropy noise evaluation requires s^ given by 


0 ) 


^ [ a dA 


(D-2) 


where is the burner exit area, and the subscript u) represents Fourier trans- 
form of the corresponding primed quantities. If the cross-correlation properties 
of the temperature fluctuations are invarient across A , it may be shown that 


Sa* » S ; - S . Y 

Es Ae aa si oi si A ai 
^ cor 


(n-3) 
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where S and Y represent the spectral and coherence estimate, respectively, of 
the subscripted quantities and 1 stands for any duussy point on the burner exit 
plane. refers to the area over which the temperature fluctuations are corre- 

lated at the burner exit plane and Is given by 

J dA(d) J dA(r) a^(t) a* G + d) 

oa e (D-4) 

where r Is a vector representing any point on the burner exit plane and d Is a 
vector representing the separation distance between any two r*. It is obvious 
from the above expressions that the area averaged quantities can be derived from 
the single point thermocouple measurements, once the correlation area as a function 
of frequency is evaluated through Eq,(D-4), The variation of the normalized corre- 
lation area with the frequency, as obtained by a single radial thermocouple tra- 
verse, is shown in Fig, D8. For reasons to be explained later, an average corre- 
lation area Is also computed. This average one Is computed from three correlation 
areas obtained throu^ radial as well as circumferential thermocouple cross-corre- 
lations without any end attachment and radial cross-correlations with the orifice 

3 

plate attached. The airflow rate and fuel/air ratio were kept constant at 9,1 m /min. 
and 0.019 respectively in the above three test runs. This average temperature corre- 
lation area is also shown in Fig. D8. It can be seen from Fig. D8 that the low 
frequency temperature fluctuations are correlated over a larger area compared to 
the high frequency fluctuations as is reasonable. Since the correlation area is 
a measure of the eddy size, it. can be said that an eddy representing the low fre- 
quency fluctuations is larger in size compared to that of a high frequency one. 

Coherence Estimates . So far, the spectral characteristics of the Individual 
signals have been discussed in detail. It will be expedient at this stage 
to look at the coherence between the various signals. As already mentioned in 
Appendix B, coherence estimates between two signals will bring out the true lin- 
ear relationship between them. Coherence analysis among the interior, near and 
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Figure D8, Frequency distribution of normalized temper- 
ature fluctuations correlation area. 





far field signals Is explained below. Various noise sources and their associated 
transmission paths are illustrated in Fig. D9. The near field signal is expected 
to contain contributions from the two interior sources, namely combustion and en- 
tropy noise sources as well as from the exterior jet noise. Hence, the coherence 
between the interior and near field signals has to be evaluated on the basis of 
a multiple input problem. In such cases, the degree of linear relationship be- 
tween any single cause and the effect is best revealed through the partial co- 
herence function analysis which cancels out the effects of the other extraneous 
inputs. The theory of partial coherence functions and the procedures for evaluat- 
ing them are described in detail in Ref. Dl. 

The results of the coherence estimates between the interior and the near 
field signals are shownlnFlgs. DIO » D13, It is important to note 
that the interior sources, namely, the interior pressure transducer and the 
thermocouple signals, may or may not be correlated. In either case, the coherence 
analysis takes care of the situation except that in the uncorrelated case, the 
expressions become simpler. Hie significance of the experimental coherence re- 
sults can be best explained with the help of some analytical expressions for co- 
herence estimates based on the model shown in Fig, D9. For the purpose of ex- 
planation, a model of an uncorrelated case will be assumed because of its simpler 
analytical expressions, although the same conclusions can be arrived at, even in 
the case of correlated inputs. For the uncorrelated interior sources, the model 
shown in Fig, D9, yields the following expressions 

Ordinary Coherence Estimate . 


Interior and near field microphones = y 
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Figure D9, Noise sources model for coherence estlmaues. 
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Hiermocouple and near field microphone <■ V 2 y 


1 



(D-5) 


Partial Coherence Estimate . 

Interior and near field microphones 
with entropy noise effect removed 



Thermocouple and near field microphone 
with interior microphone effect removed 




where S denotes the many sample averaged power spectral densities of the sub- 
scripted source. Returning to the coherence results, it can be seen from Figs. 

DIO and Dll that the ordinary coherence function level between the interior and 
near field microphones generally d''creases with an increase in the exit Mach 
number, while that between the thermocouple and the near field increases. This 
result, interpreted through the Eqs. (D-5) suggests that the entropy noise con- 
tribution available to the near field microphone starts dominating the combustion 
noise contribution to the near field microphone signal, with the increase in the 
exit Mach number. Stated otherwise, there is a gradual transition from a combus- 
tica noise dominated situation at low Mach numbers to an entropy noise dominated 
one at high Mach numbers. 

The partial coherence estimates, a more reliable procedure for the multiple- 
input problem, also confirm the above results in Figs, D12 and D13, It can be 
seen from these figures that with the increase in the exit Mach number, the par- 
tial coherenc."' estimate obtained between the thermocouple and che near field sig- 
nals, with the interior microphone effects removed, rises above the partial co- 
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D12, Partial cohsrance sstimattis betuaen lnte» 
rior and near field signals for a fixed 
airflow rate and fuel/air ratio with vary- 
ing exit terminations. 




67 




-r* 




PARTIAL COHERENCE 

INTERIOR a NEAR FIELD MICROPHONES 
( WITH ENTROPY NOISE EFFECT REMOVED ) 

THERMOCOUPLE a NEAR FIELD MICROPHONE 
( INTERIOR MICROPHONE EFFECT REMOVED) 

7.0 0.016- O.C 


1 9« 8 — 0« 0l6 — I • 0 


0.6- \ A 

\ f \ 

\ ' ^ \ 

0. 4 - ^ \ 

/\ \ / 

/ \ \ / 

0 * 2 -/ \ \ 

1 I 1 ..j. 

lOv- 500 

FREQUENCY (H::) 


Figure D13. Partial coherence estimates between inte 
rior and near field signals for a fixed 
exit termination and fuel/air ratio with 
increasing airflow rates. 





hereuce levels evaluated between the Interior and the near field microphone 
signals, with the effect of the entropy noise ranoved. This partial coherence 
result, with the help of the Eqs. (D-6) reestablishes die previous result that 
the entropy noise starts dcMninating the near field signal with the increase in 
the exit Mach number. In other ifords, the partial coherence analysis also re- 
affirms die results of the ordinary coherence function analysis that a combus- 
tion noise domination of the near field signal at low Mach numbers yields to 
an entropy noise domination at hi^ Mach numbers. It is important to note that 
in all the above partial coherence calculations the single point thermocouple 
measurements were converted into area averaged quantities using Eqs. D-3 and D-4 
which involve the ratio ^ at Fig. D8 reveals a significant 

difference in the correlation area values ccmiputed in two different ways. The 
coherence results, discussed above, are based on the correlation area obtained 
from a single radial traverse of the thermoco' c at the burner exit plane with- 
out any nozzle or orifice plate attachment. It is Informative to see how sensi- 
tive the coherence estimates are to die correl9tion area values. Referring to 
Fig. D9, the partial coherence estimate between 1 and y, wldi the effects of 
2 subtracted out by a linear least-square prediction, is given as follows, 


^ly,2 “ ^ly 


^ cor *22 ^Iv 

(l- Al) (* - ^ ^2y) 

' cor / \ cor ' 


(D-7) 


The symbols used above have been explained in the previous equations and the 
ratio appears in the above equation to obtain die area averaged quan- 

tities from single point thermocouple measurements, as explained in the Eqs. 

(D-3 and 4). For a given test run, the spectral and ordinary coherence estimates, 
after ensanble averaging and spectial smoothing, have very stable values. Hence, 
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the accurate estimation of die partial coherence estimate in Eq. (D-7) depends 
heavily on die confidence with which can be computed fcr a given A^. Fig. 

D8 compares the correlation area obtained in two different ways. As explained 
in die previous section, die average correlation area is computed from die three 
correlation areas obtained throu^ the radial as well as the circumferential 
thermocouple cross-correlations at the burner exit plane without any end attach- 
ment and the radial cross-correlations at die burner exit plane with the orifice 
plate attached. The otb<>.r one is obtained by a single radial traverse of the 
thermocouple across the burner exit plane without any end attachment. It can be 
seen from Fig. D8 that diere is about 2:1 change in the values when computed 
in two different ways. For a given set of spectral values, this change may drive 
the bracketed quantities in Eq. (D-7) co fluctuate between very low positive 
values and negative ones, thereby altering the partial coherence estimates sig- 
nificantly. Fig. E14 shows the partial coherence estimate based on this average 
correlation area for a few typical cases. A comparison of Figs. 13.2 and D14 with 
the help of Fig. 06 indicates that a decrease in the correlation area values 
generally increases the partial coherence estimates, with its effect being felt 
more on the coherence between the thermocouple and the near field signals than 
on die one between the interior and the near field microphones. Ihe above argu- 
ments lead to the fact that a sufficiently accurate determination of the tempera- 
ture correlation area distribution is vital for a reliable estimation of the par- 
tial coherence. This may be achieved by using a large number of thermocouples at 
the buimer exit plane so as to obtain more details about the spatial distribution 
of the temperature fluctuations. Ihe extraordinarily high coherence at about 100 Hz 
in the case of high Mach numbers is explained at the end of this section. It is im- 
portant to note that the coherence values in some cases, even after ensemble averag- 
ing, are found to be in the range of the statistical errors associated with the co- 
herence estimates. With a view to increasing the reliability of the results, a spectral 
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Figure D14, Partial coherence e.timates between inte- 
rior and near field signals for different 
exit terminations, using average tempera- 
ture correlation area. 
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smoothing technique has been employed In which die spectral estimates corres- 
ponding to 9 adjacent frequency bands have been averaged. This process reduces 
the frequency resolution but Improves die statistical stability. 

Ihe results of die ordinary coherence estimates between the near and far 
field oiicrophones are presented in Fig. D13. A good coherence level exists 
at low exit Hach numbers. This is attributable to the fact that at these low 
Mach numbers, the combustion noise contribution to the exterior radiated sound 
is very hi^, and, further, its low frequency nature yields to a monopole radi- 
ation pattern. However, at high exit Mach numbers, the coherence level starts 
decreasing as shown in Fig. DI5. This may be due to the fact that with the in- 
creasing exit Mach number, the jet noise starts contaminating the core noise 
radiation. The far field microdione signal is contt>.ninated to a greater degree 
than that of near field as the former receives the effect of whole jet. However, 
the near field micro{hone receives greater contribution from core noise radia- 
tion with a lesser degree of contamination by jet noise. This leads to low co- 
herence levels between near and far field signals at high exit Mach numbers. 

It can be seen from the coherence results of Fig. ElO to D15 that in the 
high Mach number cases, the coherence suddenly jumps to a high value in the 
vicinity of 100 Hz. The interior and the near field microphone spectra as well 
as the thermocouple spectra exhibit a similar btdiavior at the corresponding fre- 
quencies, as shown in Figs. D2,D3 and D5. Many possibilities for the occurance 
of this low frequency peak have been investigated. A speculation about the con- 
tribution by the burner stand or the probe to this peak has to be discounted, 
because this occurs only in the high Mach number cases and, moreover, the corres- 
ponding cold flow tests do not exhibit this trend. A suspicion on a longitudinal 
resonance phenomenon is ruled out because the resonant frequency of the combustor, 
calculated from the speed of sound within the burner rig and the length of the 
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Figure D15, Ordinary coherence estimates between near 
and far field microphones. 
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combustor, falls in the range of 400-800 Hz, depending on the nature of the 
exit termination, as supported experimentally by Fig, D2. A strong possibility 
is that there may be an acoustic coupling between the interior pressure waves 
and the fuel/air ratio. Pressure waves traveling towards the head end of the 
burner will bring out a change in the density of the air which brings a change 
in the fuel/air ratio. This changes the temperature of the fluid traveling 
towards the nozzle. Uie temperature change encountering the nozzle in turn pro- 
duces a new set of pressure waves, augumenting the old set. It is believed that 
this acoustic feed-back is mainly responsible for the observed low frequency 
(about 100 Hz) peak, since the frequency of the phenomenon should be governed 
by the flow speed and the burner length. Calculations show that this frequency 
falls in the range of 100-200 Hz, This peak, while interesting, is not too sig- 
nificant as far as the overall sound output is concerned in these experiments. 

Finally, note from Eqs, (D-6) that the partial coherence function should be 
unity if there are no intervening noise sources other than the two under con- 
sideration (S =0), A partial coherence of unity was not achieved with either 
zz 

partial coherence function for either area correlation curve. Moreover, this 
was not even achieved for runs where the jet noise contamination was low. One 
could possibly conclude, therefore, that there are some unconsidered noise 
sources which are important to the core noise problem. Unfortunately, the par- 
tial coherence function is relatively sensitive to the correlation area meas- 
urements so that this conclusion appears weak. On the other hand, the partial 
coherence function between the interior and near field microphones should not 
suffer from errors in the correlation area, and at high exit Mach number there 
is indeed low partial coherence between these two microphones. There is the sus- 
picion, therefore, that at least one other noise source is present which was not 
investigated in this program. The most likely candidate appears to be vorticity- 


nozzle interaction noise, which deserves investigation in a future program. 
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Appendix E 


The Relation Between Direct Combustion and Entropy Noise 
In Appendix C it was shown that in the absence of any containment to an 
interior microphone that the pressure transform is of the form 


= f(u)) 


dV Q cos k X 
tt) oo 


(El) 


where f(ou) is a frequency dependent transfer function, V is the combustor 
volume, k is the plane wave mode wave number and Q is the transform of 
the dimensionless heat release fluctuation. Also, it was shown that the en- 
tropy follows the equation 

52. = Xri ^ 

Dt Y 

which upon linearization and application of the Fourier transform becomes 


i(UO + v,va+v ,7a 
to -* uu -*u) 


izl 


U) 


(E2) 


Assuming the mean flow is nearly one dimensional, Eq, (E2) becomes 

^ d CT v-i • 

icua + n - — ^ + u -3 — = ^ — Q 
U) uj dx uu 


Taking a cross section average by the operation — J 

da . - IIP 

- u) , . A , d g Y-1 1 r i ,0 

dx UJ (Uqq dx Y S J u) 


(E3) 


For purpose of analysis it will be assumed that the convective term u^^da/dx 

^ • 

is negligible compared with the fluctuations in due to Q^. This appears 

justified if u is only the acoustic velocity, but theoretically the tur- 
U>oo 

bulent fluctuation (vortical part) should also be considered. More sophisti- 
cated analysis can reexamine this issue in future work. Equation (E3) may then 
be written 


d5 

dx 


CT 

u U) 
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yU S 


Q dS 
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(E4) 
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If the boundary condition is applied that there is no combustion, and, hence, 
no fluctuation in a at x = 0, the solution to Eq, (E4) may immediately be 
written as (evaluated at the exit plane) 


(U. 


YS 


P dx ^ iuif 

J -r r e qJ u I i 


a j 


dx' ! Q dS 
S 


which may be computed (since S = V in the non-dimenionalization scheme used) 



dx 

T 


1 

V 


r 




Equation (E5) may be written as 


a = g(ou) 



dv 


(E5) 


(E6) 


where g(u)) is a deterministic transfer function. It can already be seen that 
there is a remarkable similarity between Eq, (E6) and Eq, (El), This is espe- 
cially true if u)/vi <<1 in which case cos k x and the exponential terms 

in Eqs, (El) and (E6) are unity. With the nondimensionalization scheme used 
here Ui ^ n,a v ^ ^ corresponds to the physical frequency being of the same order 
of magnitude as the inverse of a fluid particle stay time In the combustor. 

In the opposite limit < < 1 there is extremely rapid phase rotation 

of the exponential factor in Eq, (E6) as the x-integration is carried out 
which would destroy any phase coherence between Eq, (El) and (E6). 

The expected situation is .consequent. y, near perfect coherence between 
p and a [except for the factor u in the denominator of Eq. (E6)] at low 

U) u>g 

frequency, changing to zero coherence at high frequency. The transition 
whould take place where oo « 1" the current experimental hardware the 

physical frequency is roughly 200 Hz, 
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5 is the cross section average transform of the entropy fluctuation, A 
'“e 

point measurement will be expected to yield lower coherence, as discussed in 
Appendix D, Tlie coherence function, for a nozzle-off case, between a point 
temperature measurement and the wall pressure fluctuation is shown in Fig, El, 
The drop in coherence at very low frequency is believed due to the hydrodynamic 
noise contamination of the interior microphone, as discussed earlier. The drop- 
off above 150-200 Hz is the transition from coherence to incoherence between 
p and a . 


U) 




This phenomenon is shown also in Fig, E2 where the cross correlation co- 
efficient between p^ and is shown. Although detailed interpretation of this 
curve is impossible because Eqs, (E6) and (El) contain complicated phase in- 
formation, it is readily observed that a) a strong cross correlation exists 
and b) the "frequency" of the cross-correlation oscillation is roughly 130 Hz, 
which is of the expected magnitude. 

The major conclusion is that above about 200 Hz it is permissible to assume 
statistical independence of combustion noise and entropy noise, whereas at lower 
frequencies they are intimately connected. This transition frequency, of course, 
depends upon the flow speed and burner size. 
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Thermocouple Time Constant 
Measurement by Cross Power Spectra 

Warren C. Sirahle* and M. Muihukrishnant 
Georgia institute of Technology. Atlanta, Ga. 



Iniroduclion 

T he measuremeni of fluctuating temperatures down- 
stream of the combustor in turbopropulsion systems is 
required to quantify the importance oi entrops noise 
generation in these systems. ' ' It is well known, howe'er, it’ ii 
most thermoi'ouples suitable for use in such a hostile cn 
vironment have response times considciably longer tli.in 
required for flat response in the audible frcquencs .angc 
Consequently, they must be compensated. Central to the com 
pensation problem is the problem of measuring the response 
lime of a g'.'en thermocouple, because it must be measured in 
the environment which it will see in use. This is ‘o because the 
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TECHNICAL NOTES 


response lime depends upon the convective and radiative en- 
vironment in which the tbermocoupie is located.'* Methods 
previously developed for measurement of the time constant 
have usually relied on measurement of the impulse response 
by quickly inserting the thermocouple into the environment, 
switching a gas stream from cold to hot, or by turning off an 
initial current through the wire and watching the temperature 
decay to the temperature of the environment.' Mechanical 
methods are adequate for time constants of the order of 100 
nisec or longer. However, for “fast" response thermocouples 
( < 10 msec), the mechanical switching methods are not 
usually fast enough to provide a true step input temperature 
change. Electrical methods, since they involve heating of the 
wire above the environment tempetature. are not satisfactory 
if the thermocouple is operating near its limit of survivability. 
The purpose of this note is to show a new method of ther- 
mocouple lime constant measurement, which at the same lime 
is nonintrusive upon the thermocouple. 

Aoalysb 

The fundamental assumptions are the equation for the AC 
output of the thermocouple, y, which is proportional to the 
apparent temperature fluctuation about the mean tem- 
perature, is eiven by y=fyz, where p is the portion of the 
signal truly proportional to the temperature flucturations. 
and c is the background electronic noise, p obeys the linear 
differential equation 

dpi, 

(I) 

at T 

Here r is the time constant, x is the true temperature fluc- 
tuation. and r is time, ij is a factor, usually very near unity, 
which may be calculated from known radiation and con- 
duction correction factors. * it is the ratio of the equilibrium 
thermocouple reading to the reading it would give if it were at 
the true temperature. The method of measurement of r is the 
issue. In the method outlined here, two thermocouples, 
denoted by I and 2. of di.'fering time constants are required to 
be placed in the flow as closely as possible to one another so as 
to sec the same temperature history. Any separation of the 
two thermocouples will yield some error, the exact nature of 
which is discussed later. 

The method assumes the availability of equipment which 
will readily process the signals by Fourier analysis. Taking the 
firaie Fourier transform of Eq. ( 1 1. 


and by definition 


T( / + loir) =j| A' 


ifA 

Y=Y + Z=~ hZ 

/ +IWT 


Here capital letters stand for 'he Fourier transform of the 
small letter quantities. The cross power spectrum and the auto 
power spectra are^ 


Si: = y, Y: 


S,i = Y, Y-, 


s., = y,r. 


which yield from Eq. (2) 


S,,.( 

' / +/WT, ' ' / -/u>T, ' 

' / + IUT/ ' ' / -tJir ! ' 

{Here, because they will cancel oiii la'er. proporlionalily eonslarn in 
ihe speciral funclionc incoiving ihc camplmp time have been omillcd 


where a * denotes the complex conjugate. If the spurious noise 
is incoherent with the signal then an ensemble average, 
denoted by a bar superscript, will yield 

/ +ur r, r, + ice! r/ - Ty ) 

+Z/T; (4) 

provided that z, is incoherent with Z:- In what follows it is 
presumed that the spurious noise for signal I is sufficiently 
small that Z,Z|< may be neglected in comparison with the 
first term in New constructing the ratio of S,, to S/>. 




/ + ul ■^T/ T_» +iw(r/ — T. ) 

l+W-T,‘ 


Consider then 


lim /? = /?, = r^/T, 


which directly gives the time constant ratio. More im- 
portantly, cc'isider the imaginary part of R, R,. 

1+Ql‘Ti' 

I^ffeientiating, and finding an extremum in this function. 
ri-f! 2u)-(T,-T..)r, • 

f+toljT, • U+U-T,')' 

where the e subscript denotes an extremum value of w. For 
T;^T., which is the reason for using two different ther- 
mocouples, there is a single extremum at 

w,. = I It, (7) 

Thus, T, may be found directly if a plot of /?, is .available. The 
primary restrictions here are that both thermocouples ^re 
seeing the same x and that the noise aoes not mask the signal 
in the vicinity ofio... 

The quality of the data can be assured by constructing the 
coherence function. This is 

, S,;S i; I 


Z,Z,’U + {-^T,) '] ^ Z,Z,‘|/-t-(a,r, ) -’l 

n'xx* r’AA* 


Z,Z',Z;Z':lt,)T,)-] [ / (wT. ) ■’) 

i/'AA* 

All the terms in b are positive, real qualities. .^ sufficient 
condition, therefore for Eq. (7) to be valid is one of > ■’ being 
close to unity in the vicinity of u , . 7 nis condition furthermore 
assures that X at position I is ihe same as \ at position 2 since 
the effect of any difference i» to inscii a spurious noise into 
one signal that is not seen at the other This coherence con- 
dition of unity is not a necessary one, howevci. All .hat arc 
necessary are the conditions leading to Eq. (5), in ihc vicinity 
of ui , . The rc.?son the coherence condition is introduced is that 
there is commercially available equipment that will readily 
compute ■> • . 

Finally. Eq. (6) is not useful as Eq. (7) because it requires a 
large value of frequency. Si-ice both signals will be falling off 
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marKedly at high frequency there is much more danger of 
dropping into the bcck^.ound nc c at high frequency as com- 
pared with u),. 


Experiment 

To illustrate the procedure two commercial Chromel- 
Alumel thermocouples of wire diameters 0.00! in. (TCI) and 
0.003 in. (TC2) were mounted less than 1/8 in. apart in the 
exhaust of a gas turbine combustor. Fixed bandwidth spectra 
(15.6 Hz) of the AC component of tem'.craturc are shown in 
Fig. 1. Also shown in Fig. I is the true temperature spectrum 
{il^X,Xi*) after correction, knowing the thermocouple time 
constant. Since the time constant goes roughly as D' where 
D is the wire diameter, TC2 has a poorer response, which is 
also shown in Fig. I . In this example the ratio of the time con- 
stant, Tj/T,. should be roughly 3 ' ' = 5.2. 

The time constant for TCI was expected to lie between I 
and 10 msec, so Finer bandwidth filtering (3.1 Hz), and z 
narrower frequency range (0-400 Hz) was chosen for time 
constant analysis in Fig. 2. In Fig. 2a, the coherence function 
is seen to be adequate within the band 10-400 Hz. in Fig. 2b 
are shown R, and R,. R, shows ih«. expetted behavior of 
monotonically going from unity to t ./r, as goes from zero 
to values much greater than I 'r.. R, shows the expected 
behavior of a maximum at 40 Hz, corresponding to 
T, = l/(2»-40) = 4msec. 

Some problems .vith this method are a) the AC component 
of temperature is usually substantially lower than tne DC 
component so that hig:< systems gains are needed and higher 
than desired background noise is usually evident, especially 
for the larger thermocouple, and b) the maximum is 
somewhat broad in R, so that, in this example, about 25S 
error in r may be expected. Problem a) is minimized by using 
thermocouples of nearly equisalent (hu' not equal) time con- 
.stants. However, it appears that problem b) must be accepted 
as a limitation of the method. 
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Appendix G 


Liner Hole Size Variation Investigation 

The special manufacture of different liners was described in Appendix A* 
These were tested in the nozzle-off configuration to determine spectral shifts 
in the condnistioc noise due to liner hole size modifications. Only the near 
field microphone was used, it being the one to sense primarily combustion 
noise. 

The results of the i^ar field spectra for several runs are shown in Figs. 
Gl - G3. The raw spectra show little, if any, change in frequency content 
with a change in liner hole size. To put this on a more precise basis, the 

3 

spectra for the 14.2 c /min air flow, 0.02 fuel/air ratio cases were inte- 
grated over the range 0-2000 Hz. The results are presented in Fig. G4; from 
this plot the "half-power" frequency may be seen, above and below which half 
of the cumulative noise exists. There is virtually no shift la frequency with 
liner hole sl^e. If anything, there is a mild shift to higher frequency with 
an increase in hole size, contrary to physical expectation. 

To put the results on an even more fundamental basis, the technique of 
Appendix C is applied to extract the h^at release spectrum for a set of three 
liners at fixed fuel/air and mass flow. Shown in Figs. G5 - G7 are the heat 
release, hydrodynamic noise, and interior microphone spectra as well as the 
Integrated (or cumulative) spectrum of the heat release function. The heat 
release spectrum is virtually unchanged with a change in liner hole size. 

The results were surprising, since a previous combustion noise study'' ’ 
showed a correladon of frequency content with turbulence properties, in- 
cluding turbulence length scale. The hole size variation was specifically 

(G2) 

carried out to modify turbulence length stale. Previous results' ' have 
shown that virtually no flow, heat release, or conbustion rate property 
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Figure G3, Near field spectra for three different 
liners. 
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affects the frequency content. This is an untenable situation. In the ex- 
periments of Ref. (G2), however ) the fuel reactivity was only varied over 
a factor of about 2:1, and the more sophisticated analysis used here was 
not applied to those results. It is believed a more corrected effort with 
a wide range of fuel reactivity should be carried out to determine if fuel 
reactivity is the determinant of frequency content. 
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